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Abstract 

In this work we consider a four compartment Brusselator system. The reaction terms of this system 
are of non constant sign, thus components of the solution are not bounded apriori, and functional means to 
derive apriori bounds will fail. We prove global existence of solutions, via construction of an appropriate 
lyapunov functional. Furthermore due to the sign changing nonlinearities, the asymptotic sign condition 
is also not satisfied, causing further difficulties in proving the existence of a global attractor. These 
difficulties are circumvented via the use of the lyapunov functional constructed along with the use of 
the uniform gronwall Lemma. We are able to prove the existence of an (L 2 (Q), H 2 (Q)) attractor for the 
system, improving previous results in the literature from [23] . The Hausdorff and fractal dimensions of the 
attractor are also shown to be finite. In particular we derive a lower bound on the Hausdorff dimension 
of the global attractor. We use numerical simulations, as well as numerical attractor reconstruction 
methods via non linear time series analysis, to validate our results. 
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1 Introduction 

The object of the current manuscript is to investigate global existence and the long time dynamics, of the 
following four compartment Brusselator type reaction diffusion system. 


du 

W7 “ aAu 
at 

dv 

m - bAv 

dw 

WT - cAw 
at 

dz 

m~ dA 2 


a — (/3 + 1) u + u 2 v + D\ (w — u ), 
/3u — u 2 v + D 2 (z — v ), 
a — (/3 + 1) w + w 2 z + D% (u — w ), 
f3w — w 2 z + D 4 (v — z ), 


(1) 

( 2 ) 

(3) 

(4) 
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in M + x $7, with the homogeneous Dirichlet boundary condition 


u = v = w = z = 0 in M + x dTl. 


( 5 ) 


We also impose suitable initial data 

u(0,x) = uo(x), v(0,x) = vo(x), w(0, x) = wo(x), z(0, x) m zo(x) in ft. (6) 

Here is an open bounded domain of class C 1 in R N , with boundary dQ . The constants a , b, c, d, a, /3 and 
Di, i = 1,2,3 and 4 are positive. The initial data are assumed to be nonnegative. The reaction terms denoted 
respectively by /, g , h and k are continuously differentiable functions on satisfying /( 0 , v , re, z) > 0 , 
g(u , 0 , re, z) > 0 , /i(iq r>, 0 , z) > 0 and k(u, v, re, 0 ) > 0 for all u, v, re, z > 0 which imply, via invariant region 
methods [16], the positivity of the solution on its interval of existence. Several authors as [ 6 ], [7], [5] , [ 12 ] 
and [16] established global existence for solutions of m-components systems (m > 2 ) with the boundary 
conditions (1.5). 

To the best of our knowledge the most recent work [23], that considers the system (1) - (5), does so under 
the following constraints 


a = c, b = d, D { = £> 3 , D 2 = D A . (7) 

This constraint is required, as the techniques in [23], require addition of the equations, to derive apriori 
bounds on the sum of the resultants. The addition causes the problematic nonlinearities to cancel because 
of (7), which facilitates the apriori estimates. In [23] Global existence, existence of a (L 2 (fl)(fl), Hq (11)) 
global attractor, and existence of an exponential attractor are all demonstrated, under the assumption via 
(7). The primary contributions of the current work are 

• We remove the assumptions via (7), and still derive global existence. 

• We remove the assumptions via (7) and show the existence of a finite dimensional global attractor. 

• We improve the regularity of the attractor, and show that it is infact an (L 2 (fl), H 2 (Q)) attractor, if 
the space dimension N < 3. 

• We present a new lower bound on the Hausdorff dimension of the attractor. 

• We use nonlinear time series analysis, to numerically estimate this lower bound. We also quantify 
temporal chaos in the system. 

• We perform numerical simulations to elucidate the chaotic dynamics on the attractor, for a case where 

(7) does not hold, that is D\ / D 2 ^ ^ D^. 

Remark 1.1 We would like to point out that the case with D\ ^ D 2 ^ ^ D 4 , as we consider, can 

lead to very interesting dynamical behavior, including aperiodic dynamics or a chaotic attractor. This is 
validated via numerical simulations in section 8. For the numerics we use Neumann boundary conditions. 
Note however, the existence results hold under these boundary conditions as well. 

We first present a brief reveiw of the relevant literature. The Brusselator, in its original form, is a system 
of 2 ODE’s that model cubic autocatalytic chemical reactions [14]. The diffusive Brusselator system is given 
by 


du 

at “ d,A " 

^ 

ot 


a — (b + 1 ) u + u 2 v , 
bu — u 2 v. 
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( 8 ) 

(9) 

( 10 ) 


This system exhibits rich dynamics, including oscillations, spatiotemporal chaos and turing instabilities 
[15, 17]. There has been a recent interest in the global dynamics of such systems, and a number of works 
have appeared to this end [21, 22 , 23]. Equations of this form pose various challenges from a mathematical 
point of view. Note, the reaction terms in (l)-(5) do not have a constant sign, so neither component is a 
priori bounded, or at least bounded in some L p - space, in order to apply the well known regularizing effect. 

From the point of view of long time dynamics, an inherent difficulty in systems of the type considered, is 
that the asymptotic sign condition in vector version 

limsupFfs) s < C ( 11 ) 

| s | —yoo 

(where C is a positive constant and F is the nonlinear term, representing the reaction), is not satisfied. 
This again is primarily due to the opposite signed terms u 2 v,w 2 z and — u 2 v,—w 2 z in equations (1) - (5) 
Usually this condition plays a key role in the dissipation process, and thus if it is satisfied, often leads to the 
existence of a global attractor, for the system. This opposite signed nonlinearity however, best represents 
the chemical process at work. From a mathematical point of view, this opposite signed coupling, causes 
extensive problems in proving existence of a global attractor. Essentially, showing asymptotic compactness 
of the semigroup in question, is not straightforward anymore. This difficulty is circumvented by making 
various apriori estimates, where the key tool used is the uniform Gronwall lemma. 

The organisation of the current manuscript is as follows. In section 2 we introduce various preliminaries 
that are required throught the manuscript. In section 3 the global existence of strong solutions is proved 
via Theorem 3.2 and Proposition 3.4. In section 4 we prove the global existence of weak solutions, and 
construct absorbing sets, in the phase space L 2 (U) via Proposition 4.4 and Lemma 4.5. In section 5 we make 
further estimates as required for the existence and regularity of the global attractor via Lemmas 5.5, 5.7. In 
section 6 we show the existence of a global attractor via Theorem 6.5. Section 7 contains results on the finite 
dimensionality of the Hausdorff and fractal dimensions of the attractor, by providing upper bounds for them 
via Theorem 7.4. Furthermore, we also present a lower bound on the Hausdorff dimension of the attractor 
via Theorem 7.5. Section 8 is devoted to numerical simulations. Here we essentially elucidate the chaotic 
dynamics of the attractor, for the case that D\ ^ ^ ^ D±. Lastly, in section 9 we use nonlinear time 

series analysis to estimate lower bounds on the attractor, as well as show temporal chaos. 

Also in all estimates made hence forth, C,Ci,i = 1,2,3.. are generic constants, and can change in value 
from line to line, and sometimes within the same line, if so required. 


2 Preliminary observations 

The usual norms in spaces ZT(U), L°°(ft) and C(Q) are denoted respectively by 

ll< = [ \n(x)fdx, (12) 

IMIoo = maxM;E)|. (13) 

oc G “ “ 

It is well known that to prove global existence of solutions to (l)-(5) [4], it suffices to derive a uniform estimate 
of \\f(u,v,w,z)\\ p , \\g(u,v,w,z)\\ p , \\h(u,v,w,z)\\ p and \\k(u,v,w,z)\\ p on [0 ;T max [ for some p > N/2. Our 
aim is to construct polynomial Lyapunov functionals allowing us to obtain L p — bounds on u\v,w and z 
that lead to global existence. Since the functions /, < 7 , h and k are continuously differentiable on IR\, then 
for any initial data in C(f2), it is easy to check directly their Lipschitz continuity on bounded subsets of the 
domain of a fractional power of the operator 

( —aA 
0 
0 

V o 


0 

-bA 

0 

0 


0 

0 

—cA 

0 


0 \ 
0 
0 

-dA ) 


(14) 
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Under these assumptions, the following local existence result is well known (see [13] . 


Proposition 2.1 The system (l)-(5) admits a unique, classical solution 
(u,v,w,z) on (0,T max [xfi. //T max < oo then 

iim {|| U (t, OIL + lb (t, OIL + \\w (t,OIL + lbLOIL) = +°° 

L / J- max 

where T max {||u 0 |L > IbolL > ILolL > IbolL) denotes the eventual blow-up time. 


3 Global existence of strong solution 

3.1 Results 

Put ^ = i%,An = ^ c ,A 24 = ^£,A 3i = Let 0,a and p be three 

positive constants such that 

e 2 > a 2 2 , (is) 


and 

where 


(0 2 - A\ 2 ) (cr 2 - A\z) > (Ais - ^12^23) 2 , 

a • v > r 2 


A 

V 

r 


{Q 2 — A 2 2 ) (cr 2 — A 23 ) — (A 13 — A 12 A 23) 2 
(6 2 — A 2 2 ) (cr 2 p 2 — A 24 ) — (A 14 — A 12 A 24) 2 

( 0 2 — A 2 2 ) (A 34 CT 2 — A23A24) — (A13 — A12A23) (A14 — A12A24) • 

Let us define three positive sequences {O r } r>0 , { (J q} q >o and {Pp } p >0 satisfying 

6 r 6 r -\- 2 2 


*?+l 


= 0 Z , 


and 


Pp+i 


CqO'q+2 _ 2 

2 ^ ’ 

^+1 


= p , 0 <r<q<p<n. 


(16) 

(17) 


Remark 3.1 We can enforce some or all of the above sequences to be increasing, and others decreasing, by 
choosing 

= C e 0 2r , °A+A = C a o 2q and ^ = C pP 2p , 0 <r<q<p<n, 

0 r CTq Pp 

where the constants Co, C a and C p satisfy 

Cq, C a and C p < 1 , or C p > 1 0 < r < q < p < n. 

Our first main result concerning the global existence of strong solutions, of the Brusselator evolutionary 
system, is as follows. 


Theorem 3.2 Let (u (£,.), v (£,.), w (£,.), z (£,.)) be a solution of (l)-(5) and let 


L n (t) = / H n (u(t,x) ,v(t,x) ,w(t,x) ,z(t,x))dx, 

Jn 


(18) 
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where 


n p q 

H n (u, v, w, z) = '^f2'f2c p C q C q 6 r cr q p p u r v q ~ r w p ~ q z n ~ p , (19) 

p=0q=0r=0 

with n a positive integer. 

Then the functional L n is uniformly bounded on the interval [0,T max ]. 

Corollary 3.3 All solutions of (l)-(5) with positive initial data in L°° (Cl) are in L°° (0,T max ;L n (Q,)) for 
all n > 1 . 

Proposition 3.4 All solutions of (l)-(5) with positive initial data in L°° (fl) are global 

3.2 Proofs 

For the proof of Theorem 3.2, we need some preparatory Lemmas whose proofs are similar to those in [7]. 
Let us denote H n in (19) by H n ( O r a q p p ), then we have 

Lemma 3.5 Let H n be the homogeneous polynomial defined by (19). Then the first partial derivatives of 
H n are given by 


&uHn — 'ClR-n—l { 0 r -\-\( 7 q-\-\Pp-\-\ ) , d v H n — nH n -i ( 0 r ( 7 q-\-\Pp-\-\) i 

^wH n — nH n _i (O r <JqPp -\-\) 5 d z H n — nH n —\ (Q r (jqPp) 


(4.1) 


Lemma 3.6 The second partial derivatives of H n are given by 

d u 2 H n — Tl (jl 1 ) H n —2 q-\-2Pp-\-2) 5 &uvH n — n (jl 1 ) H n —2 {6 r -\-\(J q -\-2Pp-\-2) 5 

duwH n = n in 1) H n —2 {^rpk<^q-\-iPp-\-2 ) •> &uzH n = n(ji 1) H n —2 (^r+i^g+iPp+i) ■> 

d v 2 H n = n {n 1 ) H n — 2 {6 r <j q -\-2Pp+2) 5 &vwH n = n (n 1 ) H n — 2 (d q+i Pp-\-2?) •> ( 4 . 2 ) 

&vzH n = ^ ir 1 ) H n —2 (@r&q-\-lPp-\-l ) 1 ^w 2 -^n = {jl 1 ) H n —2 {@ r CTqPp-\- 2) •> 

^wzH n — n {n 1 ) — n {n 1 ) H n —20 r o'qPpj r \, d z 2p[ n — n {n 1 ) H n — 2 ( O r c q pp ) 


Lemma 3.7 Let A be the symmetric matrix defined by 


a 11 a 12 a 13 
_ &12 U22 &23 

^13 <^23 «33 

ai4 <224 &34 

then 

a\i (<2h< 2 2 2 - Oi 2 ) det A = (PQ - P 2 ) 

where 

P = (ana 2 2 — ^12) (^11^33 “ clis) ~ (< 2 h <223 — ai 2 ai 3 ) 2 , 

Q = (<2 h <222 — CL12) (^11^44 — «f 4 ) — (ana24 — ai2<2i4) 2 ; 

R = (an<222 — CL 12) (<211(234 — &13&14) — (^11^23 — ^12^13) (<^11^24 ~ <^12<^14) • 


a i4 

a24 

a 34 

a44 


( 20 ) 
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proof 3.8 (Proof of Theorem 3.2) Differentiating L n with respect to t yields 


L' n (t) = I (d u H n + d„H n ^ + d w H n ^ + dzHnfffj dx 

(ad u H n Au + bd v H n Av + cd w H n Aw + dd z H n Az) dx 
n 

+ J ( fd u H n + gd v H n + hd w H n + kd z H n ) 

= / +J. 

Using Green’s formula and boundary conditions via (5), and applying lemma 1 we get 

/ n — 2 p q 

EEEM 0 ? (^ 9 P T) • Tu^-V-^z^-Pdx, 

q p=0q=0r=0 


\ Pp-\-l&q-\-l@r-\-l \ 
2 Pp+l^q+l^r 
^Pp+lVqOr 

dpp(JqO r J 


where {B rqp } , r = 0, q, q = 0, p, p = 0, n — 2 are £fte matrices defined by 


( app-\-2&q-\-2@r-\-2 2 " Pp-\-2&q-\-2@r-\-l 2 " Pp+^O q-\-l@r-\-l 


Brqp — 


_ - - 2 - 

~ 2 ~Pp+ 2 cr g+ 2 ^r+l bPp-\-2&q+2@r ~2~ Pp-\-2^q-\-l^r 

~2 Pp+2&q-\-l@r-\-l ~2 Pp+l&q+l@r+l ^Pp-\-2^q@r 

V 2~ Pp+2&q+l@r 2 Pp+l^q+l^r "2 Pp+l^q^r 


and T denotes the transposate vector 


T = (Vn, Vr, Vie, Vz)*. 



( 21 ) 


( 22 ) 


From the Sylvester criterion [10], each of the quadratic forms (with respect to Vn , Vn, Vie and \7z) associated 
with the matrices B rqp , r = 0,g,g = 0,p, p = 0,n — 2 is positive, if we prove the positivity of its main 
determinants A J rqp , j = 1, 2, 3 and 4. Far /ixed r = 0, q, q = 0,p, and p = 0, n — 2 ; iee see that 

A rqp = CLPp+2CTq+20 r +2 > 0 , 

and condition (15) implies 

A^qp = abp p+2 cr 2 +2 0 2 +1 (# 2 — A 2 2 ) > 0 . 

ITe can s/ioic by elementary calculation that 


A lqp = abcpl +2 a q+2 (jl +i e 2 r+i e r (0 2 - Aj 2 ) (a 2 - A|j) - (A 13 - A 12 A 23 y 


and this yields (16) A^ qp > 0. For the last determinant, we use lemma 3 to obtain 

0 92 - A u) A $ qp = abcdp 2 p+2 p 2 p+1 a A q+1 d 2 r+i e 2 r (AV - T 2 ) 

which implies, from (17), Af qp > 0 . Consequently we have I < 0 . 

Substituting the expressions of the partial derivatives given by lemma 1 in the second integral, yields 

J - n /EEEC.W 

q p=0q=0r=0 


Pp-\- 1 &q-\-l f @r-\-1 
Pp &q \ @r 


-f + g) + 


A + ^±1 h + k p p a q 9 r u r v q - r w p - ll z {n - 1) - p dx. (23) 

J PP 


Since we can choose, < 1 and < 1 , then 

7 v r Pp 1 


f @r +1 

v^r 



u 2 v < 0 , 


f Pp+1 
V Pp 



w 2 z < 0 , 


(24) 
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then we have 


Pp +1 °q +1 


/9p+1 h + k < C% (u + v + w + z + 1) 
Pp 


Pp 


cr a 


L+1 


-n —1 p q 


J < C 4 / ^^^C^_ 1 C«C; W r v 9 - r «; p - 9 z (n - 1) - p (w + v + w + ^ + l) 

a „——n 


dx. 


(25) 


,p=0q=0r=0 


To prove that the functional L n is uniformly bounded on the interval [0,T max ] ; first we write 

n— 1 p q 

^J2T C n-i C P C g urvq ~ r ' wP ~ 9z{n ~ 1) ~ P (U + V + W + Z + 1 ) 
p=0q=0r=0 

= L n (iX, V, W, Z) + S'n-l (^, W, *) , 


where R n (u, v , ic, 2 :) and S n - 1 (ix, v , w, z)are two homogeneous polynomials of degrees n and n—1 , respectively. 

By application of Holder’s inequality to the integrals (2b), one gets the functional L n satisfies the differ¬ 
ential inequality 

L' n (t)<C 6 L(t) + C 7 L^(t), 

j_ 

which for Z = Ln can 6e written as 

nZ ' < C5Z + C7. 

simple integration gives the uniform bound of the functional L n on the interval [0, T max ] ; t/izs ends the 
proof of the theorem. 


proof 3.9 (Proof of corollary 3.3) The proof of this corollary is an immediate consequence of theoreml, 
the preliminary observations and the inequality 


/ 


(u + v + w + z) n dx < C$L n ( t ) on [0, T* [, 


for all n > 1, where Cg is a positive constant. 


proof 3.10 (Proof of proposition 3.4) Since the degree of the polynomials in the reaction terms are three, 
then from corollary 1, there exists a positive constants C±q such that 


/■ 


(■U + V + W + Z + l) 3n dx < Cio on [0, T max [, 


for all n > 1 .The reaction terms are in L°° (0, T*\L n ($7)) for some n > —.Then from the preliminary 
observations the solution is global. 


4 Global existence of weak solutions and construction of the ab¬ 
sorbing set 


4.1 The (v — z) L p absorbing set 

We first recall the following definition, 
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Definition 4.1 (Bounded absorbing set) A bounded set B in a reflexive Banach space H is called a 
bounded absorbing set if for each bounded subset U of H, there is a time T = T(U), such that S(t)U C B for 
all t > T. The number T = T(U) is referred to as the compactification time for S(t)U. This is essentially 
the time after which the semigroup compactifles. 


Consider, for a positive real p > 2 , the following functional 

K p (t)= [ ( v p + 5z p )dx , 
do 

where p > 2 is a positive real. By differentiating K p with respect to £, we get 


K' p (t) = -p{p- 1 ) J (av p ~ 2 |Vv | 2 + c5z p ~ 2 |Vz| 2 ^ 


dx 


f [v p 1 (flu — u 2 v — D 2 v +D2z) j t5z p 1 [/3w — w 2 z + D^v — D^z)] dx 
do 


+P 
= / +J. 


(26) 


where 


Using Poincare inequality, we obtain 

I<-C lP (p-l)J (v p ~ 2 \Vv\ 2 +6z p ~ 2 \Vz\ 2 flx-2C 2 J ( v p + 5z p )dx , 

C\ = min (a, c), C 2 = Ai min (a, c ) 

J = p f [- v p ~ 2 (fiuv — u 2 v 2 ) + Sz p ~ 2 (flwz — w 2 z 2 ) — (v — z) (D 2 ^ p_1 — cD 4 2 ) p_1 )] dx 
do 

By choosing S = we get 


J = p 


I 

do 


0* 


p-2 


x 7 pm; — U V‘ 


! ) + 2 (fiwz — re 2 z 2 ) — D 2 (v — z) 2 ^^v q z p 2 q 

q =0 


dx 


— PPf [ \ vP 2 + 2 ] dx. 

do 

Applying Young inequality, we get 


J < C 3 + C 2 [ [v p + dz p ] dx. 
do 


Finally, we have 


J [v p + Sz p ] dx + Cip(p - 1 ) J (v p ~ 2 \Vv\ 2 + 6z p ~ 2 \\7z\ 2 flx 


d_ 

dt _ 

+C 2 [ [v p + 5z p ]dx<C 3 , 
do 


which gives 


[ [v* + 6z p ] 
do 


dx < e ° 2t f [vq + dzg] 
do 


, ^3 

<fa+ ci- 


(27) 


(28) 
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Remark 4.2 Note that the (v, z) component of the solution to the Brusselator system gets into an L p ball, 
for data in L 2 , due to the almost instantaneous regularisation of the system, we have that the weak solution 
becomes a strong solution that belongs to Hq, and hence to L p ,p < 6. Thus the data without loss of generality 
can be taken in Hq . 

4.2 The L 2 absorbing set 

In this subsection, we will prove the existence of an L 2 absorbing set for the semi group of the Brusselator. 
We use the functional L n given in section 3, for n = 2 : 

1/2 (t) = / H 2 (u (t, x) , v (t, x) , w (t, x) , z (t, x)) dx, (29) 

Jn 

where H 2 is is given by (19) for n = 2 


H 2 {u, , v, w, z) = OqctoPoz 2 + 20 o a () p 1 wz + 20 o a 1 p 1 vz + 20 1 a 1 p 1 uz + 0 o aop 2 w 2 , , 

+ 2 O 0 O- 1 P 2 VW + 20iCFip2UW + 0 O (T2p2V 2 + 20 1 a 2 p2UV + O 20 - 2 P 2 U 2 ^ 

After differentiating the functional L 2 with respect to the time t and following the same steps as in section 
3 taking into account the positivity of the matrix £>000 given by (22) and which is in this case a constant 
matrix, we get 


L' 2 (t) + 2 C 4 J (|Vw | 2 + |Vv | 2 + |Vw | 2 + |V^| 2 ) dx 

< / R(u (t,x) ,v (t,x) ,w (t,x ), z (t,x)) dx, 

Jn 


(31) 


where C 4 is a positive constant and R is a polynomial of degree two and which the form will be determined 
later. Using Poincare inequality, we obtain 


L' 2 (t) + C 5 f n (u 2 + v 2 + w 2 + z 2 ) dx + C A f n (|Vu| 2 + |Vv| 2 + |VH 2 + |V^| 2 ) dx 
< f Q R (u,v,w, z) dx. 

From (23)-(25), the polynomial R has the following form 

R = Ri (u, w) + R 2 {v, z) + R 3 (u, v , w, z) , 

where 

Ri (u, w) = ( a\u 2 + a 2 uw + asw 2 ) , 

R 2 (v, z) = (/3 !V 2 + / 3 2 VZ + p 3 z 2 ), 
i? 3 (u, V, w, z) = (71 U + J 2 W + 73) ( 5 iV + 5 2 Z + S 3 ) 


(32) 


(33) 


(34) 


and where the coefficients are given in (24). By application of Young inequality, we can find a positive 
constants Cq and C 7 such that 


and 


f R 2 {v, z)dx <C 6 [ [ v 2 + 5z 2 ] dx , 

Jn Jn 

[ Rs (u, ,v,w, z) dx < ^ f {u 2 + w 2 ) dx + C7 f (v 2 + Sz 2 ) dx. 
Jn 2 j 0 Jn 


If we prove that 
then (32) will become 


Ri (u, w) < 0, 


(35) 

(36) 

(37) 
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( 38 ) 


L' 2 (t) + ^ f {u 2 + v 2 + w 2 + z 2 ) dx 

- Jn 

+C 4 J (|Vu| 2 + |Vv| 2 + |Vtc| 2 + |W| 2 ) dx 

< (Ce + Ct) [ (y 2 + 5z 2 ) dx, 

Jn 

then by multiplying (27), for p = 2, by 2( ° 6 ( J ^ 7 ) and adding to (38). We shall get after simplifications 


1 / 2 (t) + CsK' 2 + Cg f ( u 2 + v 2 + w 2 + z 2 ) dx 
Jn 

+C 10 j (|Vu | 2 + |Vu | 2 + |Vic | 2 + |V,z| 2 ) dx < Cu, 


where the constants are positive and independent on the initial data. 

To prove (37), we have from (23) for n = 2 

- ai = pi [a( (6>i (p + 1) - p) + aitfiDi - £> 3 ], 

—a 2 = p[ [ctq ( 0' o (/3 + 1) — /3) + (o-qOq — criOi) Di + (/3 + 1)] — /3, 

-«3 = -pi - AO + pi (/? + 1) - /?, 


where 


a-+i 

9 r 


= 9 


r & q +i _ Pp+i _ p' r = 0, g, q = 0,p, p = 0 and 1. 


07, 


The coefficient 0^3 <0, if 


/ /// < P + 1 + A 

a cro < 


Pp 


and 1 > pi > 


P 


D\ ri —<Tq0qZ?i + (/3 + 1 + D 3 ) 

The inequalities (42) are satisfied if we choose 9' 0 and a' 0 satisfying 

I+-D3 


@0 a 0 < 


D 1 


Then a\ < 0 under the following conditions 

P 


9[> 


Finally 0.2 < 0, if we choose 


P + 1 + D 1 

a' 0 < 


and cri > — 


D 3 


9[(p + l + D 1 )-p- 

P + 1 


and 


1 > pi > 


— 0q [/? + 1 + (1 — ct9) D\\ + P 

P 


K Wo [P + 1 + (1 — a @) D\] — P} + (P + 1)] 

Inequalities (45) are satisfied under the following condition 

1 


<T n < 


~9'o [P + 1 + (1 — 0 $) D\ ] + ft 


(39) 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 


Since we should have < 1 and Po,pi < 1 to eliminate the terms u 2 v and w 2 z and get (24), then (43), 

(44) and (45) are satisfied if we can choose a' 0 sufficiently small and cri sufficiently large. But there is no 
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conditions on erg and a[ only the report = a 2 which should be sufficiently large to satisfy conditions (16) 
and (17). So we can choose cr' Q sufficiently small and a 2 sufficiently large to get a[ sufficiently large. Then 
(37) is satisfied for appropriate constants. 

To get the absorbing set from (39), we should remark that from the definition of the functionals L 2 and 
IT 2 , we can find a positive constants C 12 and C 13 such that 


C 12 / ( u 2 + v 2 + w 2 + z 2 ) dx 

J Q 

< L 2 (t) + C s K 2 (t) 

< C 13 / ( u 2 + v 2 +w 2 + z 2 ) dx , t > 0 . 

Jn 


(48) 


Then (39) becomes 


j t [.L 2 (t ) + C 8 K 2 (t)} + C 14 [L 2 (t) + C 8 K 2 (t )] 

+C15 J (|Vw | 2 + |Vt >| 2 + |Vw | 2 + |V^| 2 ) dx < Ci 6 , 

(49) 


where all constants are independent of the initial data. 

The inequalities (48) and (49) together give 

||(u, ,v,w,z)\\ L2(sl ) < e~ Cl7t ||(uo,,^o,w 0 , 2 o)|| L 2 ([2 ) +Ci 8 , t > 0, (50) 

where the constants C 17 and C\g are strictly positive. 

4.3 Global existence of weak solutions 

In this subsection we shall prove the global existence of the weak solutions of the Brusselator evolutionary 
equation (l)-(4) for spatial dimension N = 1, 2 and 3. We have the following standard local existence in 
time of weak solutions analogous to that given in section 2 concerning the local existence in time of strong 
solutions 

Proposition 4.3 For any given initial data (uo, Vo, rco, zq) in [L 2 (fl)] 4 ; there exists a unique, local weak so¬ 
lution (u,v,w,z) on (0,T max [xfi of the Brusselator evolutionary equations (l)-(4) with boundary conditions 
(5), which becomes a strong solution on (0,T max ). If T max < 00 then 

lim {\\u(t,.)\\ 2 4-]\v(t,.)\\ 2 + \\w(t,.)\\ 2 + \\z(t,.)\\ 2 } =+00 (51) 

I/ -I max 

where T max {||^o || 2 , 11 ^ 0 11 2 ? ||^o II 2 »11 ^ 0 11 2 } & en °tes the eventual blow-up time. 

To prove the Proposition we find estimates, similar to those presented on the two above subsections, 
on the Galerkin approximate solutions of the initial value problem (l)-( 6 ), for more details see [ 2 ] and [ 8 ]. 
Moreover, estimation (50) shows that the local solution given by the above proposition will never blow up 
in [L 2 (Q)] at any finite time and it exists globally. We have 

Proposition 4.4 All solutions of the Brusselator evolutionary equations (l)-(4) with boundary conditions 
(5) and initial data in [L 2 (Q)] are global. 
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Also, from the estimate via (50) we see there exists a time 


( n in (|«(0)|i + K0)|i + K0)|i +1^0)11)^ 

ti = max I 0 , —-- — -- I , 

such that for times t>t\ the following uniform estimate holds 

M! + Mi + Ml! + M! - ~ft ~(i + CM) < c. 

U-12 

Here C is independent of time and initial data. Thus we deduce the following result 


( 52 ) 


(53) 


Lemma 4.5 There exits a positive constant C, independent of time and initial data, such that the ball 
B 0 = fu,,v,w,z) e [L 2 (Q,)] A : ||(u, ,v,w, z)\\ L2(n) < c|, 
is an absorbing set of the semiflow associated to the Brusselator evolutionary equations. 


5 Global attractor 

In this section we aim to prove the existence of a global attractor for (l)-(5). We restrict our selves to spatial 
dimension N = 1, 2, 3. We will use the following phase spaces, 

h = L 2 (n) x L 2 (n) x L 2 (n) x L 2 (n). 

y = Hq(ci) x Hpn) x Hpn) x Hpn), 
x = H 2 (n) n Hpn) x H 2 (n) n Hpn) x h 2 (o.) n Hpn) x H 2 (n) n Hpn). 

Recall the following definitions 

Definition 5.1 Let A C H 2 (Vt). Then A is said to be a (H,X) global attractor if the following conditions 
are satisfied 

i) A is compact in X. 

ii) A is invariant, i.e, S(t)A = A, t > 0 
Hi) If B is bounded in H then 

distx(S(t)B, A) —y 0 , t oo. 

In order to prove the existence of a global attractor we are required to show: 

i) There exists a bounded absorbing set in the phase space. 

ii) The asymptotic compactness property of the semigroup in question, [20], [18]. 

These are defined next 

Definition 5.2 (asymptotic compactness) The semi-group {S(t)} t>0 associated with a dynamical sys¬ 
tem is said to be asymptotically compact in if for any {fo,n}^ = i bounded in L 2 , and a sequence of 

times {t n —>• oo} , S(t n )fo, n possesses a convergent subsequence in H 2 (Q). 

The existence of a bounded absorbing set in the phase space L 2 (Ll) has been shown via Lemma 4.5. 

Now we can integrate (49) in the time interval [£, t + 1] to yield, 
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( 54 ) 


(I V^ll + |Vv|l + IV-^II + iv^ll) 

i: 


<\u(t)\l + \v(t)\l + \ W {t)\ 


<t) I 


Cds. 


Thus for t > t\ we see that 

rt +1 


rt-\-l pt -\-1 

J \Vu\lds< j (|Vu|! + |Vv|| + \Vw\ 2 2 + \Vz\l) ds < C 


We multiply (1) by —A u and integrate by parts over Q to obtain, 
2* + “ |a “ g 

(a — (b + l)iz + + Th(u; — u))(—Au)ch 


n 


< a 


[ lA^ldx— f ((6 + l)\Vu\ 2 \dx. + f |'u 2 ||u||A'u|dx 

Jn Jn Jn 

/ \w — u\\Au\cht 
Jn 

| 2 ) + [ \u 2 \\v\\Au\chL 

Jn 


in 
T D\ 


< l\Au\l + \n\^ + ^\Au\ 2 2 + ^(\u\ 2 2 + \w 


(55) 


(56) 


We now focus on 


< 

< 

< 

< 



(57) 


The above follows via Young’s inequality with epsilon, the embedding of Hq ^ 
estimates on the L 6 norm of the v component via Lemma 5.4, we obtain, 


L 6 , and the uniform 


J t (iVizli) + ^\Au\l < C (IVtill) 2 + (Hi + Hli) (58) 

This yields 

J t (|Vu|i) < C (|V«|1) 2 + ^ (\u\l + Mi) (59) 

Now we recall the uniform gronwall lemma 


13 






Lemma 5.3 (Uniform Gronwall Lemma) Let /?, £ and h be nonnegative functions in Lj oc [ 0, oc;M). Hs- 
swme that f3 is absolutely continuous on (0, 00 ) and the following differential inequality is satisfied 


< CP + h, for t > 0. 

If there exists a finite time T > 0 and some r > 0 such that 

rT+r rT+r rT+r 

/ ((r)dr < A, / /3(r)dr < H, and / h(r)dr < C , 

JT JT JT 

for any t > T, where A , 5 and C are some positive constants, then 

/3(£) < f -j- C 1 e^, V £ > T -\- r. 


(60) 


(61) 


(62) 


Thus the use of the uniform Gronwall lemma with T = £ 1 , r = 1 and 

CD 2 

P = |V«H,C= |Vu| t,h= (| u\ 2 2 + \w\ 2 2 ) , 

and the estimates via (55), Lemma 4.5 yield the following uniform estimate 

CD 2 

\Vu\l < (C + -^)e c < C, Vt>t 2 =t 1 + 1. 

The estimates for the gradients of the v,w,z components is made similarly. We can thus state the following 
Lemma 

Lemma 5.4 Let u,v,w,z be solutions to (l)-(5) with (uo, uo, wo, zo) G L 2 (fl). There exists a time £ 2 , and 
a constant C independent of time and initial data, and depending only on the parameters in (l)-(5), such 
that for any t > £2 the following uniform a priori estimates hold: 


(63) 


(64) 


|V«|| < C, |Vu|l < C, \Vw\l < C, |Vz|i < c, 


(65) 


5.1 Integral in time a priori H 2 estimate for u 

Now we can integrate (58) in the time interval [£,£ + 1], to obtain 


nt -\-1 

/ |A u\\ds 

r t+1 f o cd 2 \ 

< |vu(t)H + J (c (|Vti|l) 2 + (Hi + Mi) J ds <c,vt> t 2 


( 66 ) 


This follows from Lemmas 4.5 and 5.4. 

The same method is used on the other components v,w,z to yield similar estimates as above. Here the 
£2 depends only on the L 2 norm of the initial data, and C is independent of time and initial data. 
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5.2 Uniform a priori H 2 estimate for u 


In this section we make higher order estimates for the components of the solutions. In these estimates 
it is assumed that u and A u satisfy the same Dirichlet boundary conditions, which is true of all solution 
components. This can be shown rigorosly [20]. We multiply Equation (1) by A 2 u and integrate by parts 
over to obtain 


1 d\ Alfc|2 | „\x~7 ( a „. m 2 


2 dt 


+ a|V(Au)|^ 


< 


(/ —V(a — (b + 1 )u + u 2 v + Di(w — u)) • V(A u)d^j 

/ ((b + l)Vu — u 2 Vv — 2 vuVu — Di(X7w — Vu)) • V(A u)dx 
Jo 

[ ((6+l)|Vw||V(Aix)|dx+ [ |^ 2 ||V^||V(Aii)|dx 
J o J o 

+2 [ H|u||Vu||V(Att)|dx + .D 1 [ |Vw||V(Au)|dx+ [ |Vu||V(Au)|ds 
J o J o J o 


(67) 


We now focus on 


< - 


< - 


< 


< 


f |M 2 ||Vv||V(A«)|dx 
Jo 

f \u 4 \\\/v\ 2 dx+^ [ |V(Ait)| 

Jo ° Jo 

2 1 

J^(u 4 )idxj ^(\Vv \ 2 ) 3 dxj +IJ q MAu)\ 

2 1, 

[ (u 6 )dx\ 3 ( f |Vv| 6 dxV [ |V(Au)| 2 dx 
Jo / \Jo ) ° Jo 

C\Av\\ + ^ f |V(Au)| 2 dx 
° JO 


2 
a 

2 

a \J o 
2 
a 


2 dx 


dx 


( 68 ) 


This follows via the use of Holder, Cauchy-Schwartz and Young’s inequality. Along with the Sobolev 
embedding of H 2 ^ W 1,6 , and the uniform L 6 estimates via (28). We next handle 

[ MM|Vti||V(ArO|dx 

Jo 

< C [ \u\ 2 \v\ 2 \Vu\ 2 dx+ | [ |V(Au)| 2 dx 
Jo 8 Jo 

c(^J (u 6 )d x y (^J H 6 ck ) 3 (y |VM| 6 dxy + | J |V(A«)| 2 dx 

C| Au l 2 + ^ [ |V(Au)| 2 dx 
8 Jo 

(69) 


< C 


This also follows via the use of Holder, Cauchy-Schwartz and Young’s inequality. Along with the Sobolev 
embedding of H 2 ^ W 1,6 , and the uniform L 6 estimates via (28). Thus we have, 
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^+*|V(A»)H 

< C\Vu\l + ^|V(A«)|1 + CilVHl + ^|V(Au)|l 

15 Id 

+ C 2 \Av\l + ||V(A«)|1 + C 3 \Au\l + ||V(Au)|l 


This follows via Cauchy-Schwartz inequality and the earlier estimates. Thus we obtain 


+ ^|V(A«)|i < C\Vu\l + Ci|Vtc | 2 + C 2 \Av\l + C 3 \Au\ 2 2 
Now similarly we can multiply Equation ( 2 ) by A 2 v and integrate by parts over Q to obtain 


1 d\Av\% 

2 dt 


■b\V(Av)\l 


(/ —\7(bu — u 2 v + D 2 (z — v)) ■ V(A v)dx^ 

< [ (6|Vu||V(Au)|dx+ [ |it 2 ||Vu||V(Au)|dx 

J CL J CL 

! [ |w||«||Vu||V(Aw)|dx + D 1 f |V^||V(Au)|dx+ f |Vu||V(Au)|d: 
JCL J CL J CL 


!CL 
+2 


In much the same manner as earlier we can derive 

+ ^|V(Au )| 2 < C 4 |Vu | 2 + C 5 |Vu | 2 + Ce\\7z\l + C 7 \Av\ 2 2 + C 8 \Au\ 2 
Now using the embedding of H 3 ^ H 2 , and adding up (71) and (73) we obtain, 


- (|Au|| + |Au| 2 ) < <7 (|Vu | 2 + |Vu | 2 + |V^H + \Vz\ 2 ) + C\ (|Au | 2 + |Au| 2 ) 
We know via ( 66 ) and lemma 5.4 that 

rt +1 rt +1 

\Au(s)\lds <C, 

It Jt 


pt -\-1 /*£+1 rt -\-1 

J \Au(s)\\ds < C, J \Av(s)\ 2 2 ds <C, J C\ds < Ci, V£ > £ 2 , 


and 


pt+% 

J (|V«(s)H + |Vw;(s)|| + |Vu(s )| 2 + \Vz(s)\ 2 )ds <C,Vt> t 2 , 


Thus application of the uniform Gronwall lemma with 


(70) 

(71) 


(72) 

(73) 

(74) 

(75) 

(76) 


P = |Au|! + |Au| 2 ,C = C 3 ,h = C 2 (|Vu | 2 + |Vu>|i + |Vu | 2 + \Vz\ 2 ),r = 1 ,T = t 2 , 


(77) 


and similar estimates for the other components, tells us there is a time £3 = £2 + 1 , such that the following 
lemma is true 
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Lemma 5.5 Consider (l)-(5). For any solutions u,v,w,z to the system, there exists a constant C indepen¬ 
dent of time and initial data, and a time £3 = £2 + 1 , such that the following estimates hold uniformly, 


\u(t) \ H 2 <C, V£ > £3 
\v(t) \h 2 <C,Mt> £3 
|w(£)| h 2 < C, V£ > £3 
\z{t)\H 2 < C, V£ > £3 


Thus the existence of a bounded absorbing set in H 2 (Q) has also been established. 

5.3 Uniform a priori estimates for || 

From (1) via brute force we obtain 


< 

< 


du 

dt 


— J [d\ Au — a — (/3 + 1) u + u 2 v + D\ (w — u )) 2 dx 

C (\Au \ 2 2 + \u\l + \v\j + \u\l + \w\ 2 2 ) 

C. 


(78) 


This follows via lemma 5.5 and the apriori L p bounds on the solutions. Note since we have uniform H 2 
estimates now via lemma 5.5, the L 8 bound on u follows via the embedding of H 2 c —>> L 8 . Similar estimates 
can be derived for the v,w,z components. We can now state the following Lemma, 


Lemma 5.6 Consider (l)-(5). For any solutions u,v,w,z of the system there exists 
pendent of time and initial data, and a time t\ such that the following estimates hold 


du 

~8t 

dv 

dt 

dw 


< C, V£ > £3 


dt 

dz 

dt 


<C, V£ > £3 


< C, V£ > £3 


< C, V£ > £3 


of a constant C inde- 
uniformly 


We next make an integral in time estimate on We take the partial derivative w.r.t £ of (1) and 
multiply the resulting equation by and integrate by parts over Q to obtain 
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1 d 

du 

2 

+ d\ 

fdVu\ 

2 dt 

dt 

2 

\dt ) 


< 


f ( . r . ^ -.dii 2 dv du dw\ f du\ 

Jn \^ +1 + Dl ^ + u ^7 + vu ~^7 +- D i^T ) ( ^7 ) dx 


dt 


< C 


< C 


< c 


< c 


< c 


du 

2 

+ 

dw 

dt 

2 

dt 

du 

2 

+ 

dw 

dt 

2 

dt 

du 

2 

+ 

dw 

dt 

2 

dt 

du 

2 

+ 

dw 

dt 

2 

dt 

du 

2 

+ 

dw 

dt 

2 

dt 


dt 

f 2 dv du f 

la" dtdt dX + Jn 


dt 


uv 


dt J V dt J 

du N 
dt 


dx 


■Cli 


■ C|Au 


dv 


du 

dt 

2 

dt 


+ C|u|oob|, 


9?/ 

dt 


dv 


du 

dt 

2 

dt 


+ C'|Aiz| 2 |Ai;|2H! 


) +c ( 

dv 

2 

du 

2 \ 

du 

dt 

+ 

2 

dt 

J +c 

dt 


+ 


du 

dt 


Thus integrating the above in the time interval [£, £ + 1 ], for £ > £3 we obtain 


< 

< 



( 79 ) 


This follows via lemma 5 . 5 , lemma 5.6 and the embedding of id 2 (Q) ^ L°°(Q). 

We will next make a uniform in time estimate for where the previous estimate will be used. We 

take the time derivative of (1), then multiply through by — and integrate by parts over Q to obtain 


< 


< 


1 d 

2 dt 


d\ 7 u 


dt 


+ d\ 


A (to 
\dt 
du 


(-(/3 + l + D 1 )—+u —+vu—+D 1 —) -A — 


/ 

Jn 

„ f , , n , —. \ dui 2 dv du _ dw . 2 d\ 

C L { ~ {13 +1+ Dl> at + "u + ™ m +D 'm >+ T 


a du 

A m 


dx 

2 


< C 


< c 


du 

2 

dw 

2 

. i i O i i O 

du 

2 

dv 

dt 

+ 

2 

dt 

+ MsMs 

2 

dt 

+ \ u 8 

4 

dt 


du 

2 

die 

2 

. i i O i i O 

d\ 7 u 

2 

i i 4 

dVu 

dt 

+ 

2 

d£ 

+ m 8 m 8 

2 

dt 

+ Ms 

2 

dt 
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This follows from the product rule for differentiation, Cauchy-Schwartz inequality and the Sobolev em¬ 
bedding of Hq(Q) L 4 ( Q). Now using Poincaire’s inequality, and the estimates via lemma 5.5, we obtain 


d 

dt 

< C 


dVu 


dt 

dVu 


dt 


C 


du 

dt 


dw 


dt 


Ci 


dVv 


dt 


We will now derive a uniform estimate for via the uniform Gronwall Lemma, [20]. 

Now via the same methods in deriving (79) we can obtain 


rt+l 

Jt 


dWv 


dt 


ds < C, V t > £3. 


Thus via the use of lemma 5.6 and application of the uniform Gronwall lemma with 




dVu(s 


dt 


C = C, r = l,T = £ 3 


and 


h = C 


( 

du(s) 

2 

+ 

dw(s) 

2 

+ Ci 

d\7v(s) 


dt 

2 

dt 

2 

dt 


we obtain 


(80) 

(81) 

(82) 


Lemma 5.7 Consider (l)-(5). For any solutions u,v,w,z to the system, there exists a constant C indepen¬ 
dent of time and initial data, and a time £4 = £3 + 1, such that the following estimates hold uniformly, 


d\7u(t) 


dt 

dVvft) 


< C, V£ > £4 


L 2 


dt 
dVw(t) 


<C, V£ > £4 


L 2 


dt 

dVz(t) 


<C, Vt> £4 


L 2 


dt 


< C, V£ > £4 


L 2 


The estimates for the v,w,z components follow similarly. 

6 Existence of global attractor 

In this section we prove the existence of a global attractor for (l)-(5) 

6.1 Preliminaries 

Recall the phase space H introduced earlier 

H = L 2 (Q) x L 2 (Q) x L 2 (tt) x L 2 (Q). 
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Also recall 

Y = Hq(Q) x Hq(Q) x Hq(Q) x Hq(Q), 
x = H 2 (n) n Hq(Q) x H 2 (n) n H^(n) x h 2 (ll) n H^(n) x h 2 (ll) n H%(n). 

Also recall that if A is an (H, H ) attractor, then in order to prove that it is a an (iL, X ) attractor it suffices 
to show the existence of a bounded absorbing set in X as well as demonstrate the asymptotic compactness 
of the semi-group in X, see [20]. We first state the following Lemma. 

Lemma 6.1 Consider the reaction diffusion system described via, (1)- (5). There exists a (H,H) global 
attractor A for this system, in space dimension N < 3, which is compact and invariant in H and attracts 
all bounded subsets of H in the H metric. 

proof 6.2 The system is well posed via proposition 3.4, hence there exists a well defined semi-group {S(t)} t>Q 
for initial data in L 2 (Q). The existence of bounded absorbing sets in H and Y follow via the estimates derived 
in Lemma 5.4 . Furthermore the compact Sobolev embedding of 

Y ^ H 

yields the asymptotic compactness of the semi-group {S(t)} t>Q in H. The existence of an ( H,H ) global 
attractor for the model now follows. 


6.2 Asymptotic compactness of the semi-group in X 

We next demonstarte the asymptotic compactness of the semigroup in X. Attempting this directly is 
quite cumbersome, as it will involve essentially trying to make apriori H 3 (Ll) estimates, and then use the 
compactness of H 3 (Ll) H 2 (Ll). We use a more elegant method, where we exploit the form of the equation. 
We show the analysis for u , the other variables follow similarly. Our strategy is to rewrite ( 1 ) as 

du 

aAu = — — (a — (/? + 1) u + u 2 v + D\ (w — u)) (83) 

We will demonstrate that every term on the right hand side of (83) converges strongly in L 2 (Ll). Thus we 
obtain that A u converges strongly in L 2 ( fi), which will imply via elliptic regularity the strong convergence 
of u in H 2 (Ll). Since the same method yields the strong convergence in H 2 (fl ), for the other components 
v,w,z , the asymptotic compactness in X follows. We state the following Lemma 


Lemma 6.3 The semi-group {S(t)} t>Q associated with the dynamical system (l)-(5) is asymptotically com¬ 
pact in X. 


proof 6.4 Let us denote u n (t ) = S(t)uo jn and U(t n ) = t=t n - We have that 

aA %(t n ) = (U (t n ) O (/5 T 1) U n ft n ) T ^ n (^n)^n(^n) H - F) i ( yj n {tn) ^n(^n)) • 

Via lemma 5.7 we have for t>t/± 


dXu 

dt 


< C. 

2 


so for n large enough t n > £4 and we obtain 


dXu n 


dt 

2 


Thus for n large enough t n > £4 and we obtain 

|VC/ n | 2 < c. 
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Also via Lemma 5.4 we have for t>t±, the estimate 

|V ^| 2 < C. 

These uniform bounds allow us to extract weakly convergent subsequences, 

U n (t n ) —>• U weakly in Hq(Q). 
u n (t n ) u weakly in Hq(Q). 

Thus from classical functional analysis theory, see [20], and the compact embedding of Hq (ft) c — L 2 (Ll), 
we obtain 

U n (t n ) U strongly in L 2 (fA), (84) 

u n (t n ) —^ u strongly in L 2 (0), (85) 

Now recall the form of the truncated reaction term, 

fn(u n , V n , W n , Z n ) = a - (/3 + 1 ) U n (t n ) + U 2 n {t n )v n (t n ) + Di (w n (t n ) - U n (t n )) 

The convergence of the linear terms in L 2 (Ut) is standard, as we have uniform Hq(Q) estimates via lemma 
5.4- In order to show convergence of the nonlinear component of f n note 

lim Wu^Vn — u 2 v \\ 2 = lim | \u^v n - u 2 v n + u 2 v n - u 2 v\| 2 

n—t oo n—t oo 

< lim \\u 2 n v n - u 2 v n \\ 2 + lim | \u 2 v n - u 2 v\| 2 

n—> oo n—»• oo 

< lim |KHl|K+«||i|K-«||l+ lim ||u||i||w„ - v\\\ 

n—>- oo n—> oo 

< lim C , ||v n || 4 ||u n + u||l||Vu n - Vu||| + lim C*||w|I 4 I|- Vv||| 

n—>• 00 n—>• 00 

< lim C||Vw„-Vu||i + C||Vv„-Vv||i 

n—)>oo 

^0 

T7&e convergence follows as we have uniform H 2 (Q) estimates via lemma 5.5, hence strong convergence of 
the components in Hq(Q). So we obtain 

fn(u n ) f(y) Strongly in L 2 (VL). ( 86 ) 

Using the convergences via (84), (85) and ( 86 ) we obtain 

A u n —>■ A u strongly in L 2 (Ul). 

However this implies via elliptic regularity that 

u n u strongly in H 2 (Ll). 

This proves the Lemma. We can now state the following result 

Theorem 6.5 Consider the coupled Brusselator system described via (l)-(5). There exists a (H,X) global 
attractor A for this system, in space dimension N < 3, which is compact and invariant in X and attracts all 
bounded subsets of H in the X metric. 

proof 6.6 The system is well posed via proposition 3.4, hence there exists a well defined semi-group {S(t)} t>Q 
for initial data in L 2 {Ct). We already have the existence of an (H,H) global attractor via lemma 6.1. The 
estimates derived via Lemma 5.7 give us the existence of bounded absorbing sets in X . Lemma 6.3 proves 
the asymptotic compactness of the semi-group {S(t)} t>Q for the dynamical system associated with (l)-(5), 
in X. Thus the theorem is proved. 
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7 Finite Dimensionality of the Global Attractor 


In this section we show that the Hausdorff and fractal dimensions of the global attractor for the reaction 
diffusion system (l)-(5), is finite. Recall 

7.1 Upper bound on the Haursdorff and fractal dimension of the global attrac¬ 
tor 

Definition 7.1 (Fractal dimension) Consider a subset X of a Banach space H. If X is compact, the 
fractal dimension of X, denoted df(X), is given by 


df( x ) = lim SU P ~ ■ (87) 

<^o log{- t ) 

Here N(X,e) denotes the minimum number of closed balls of radius e, required to cover X. Note df(X) can 
take the value +oo. 

Definition 7.2 (Hausdorff dimension) Consider a subset X of a Banach space H. If X is compact, the 
Hausdorff dimension of X, denoted dn(X), is given by 

d H (X) = inf {d : H d (X) = 0} . (88) 

Here 

H d {X) = \im/j,(X,d,e). (89) 

e —^0 


Where 


/i(X, d , e) = inf rf : r* < e and X C U iB{x^ rf) j . (90) 

and B(xi,ri) are balls with radius ri. 

We recall the following Lemma from [20], which will be useful to derive the requisite estimates. 

Lemma 7.3 If there is an integer n such that q n < 0 then the Hausdorff and fractal dimensions of A, 
denoted dn(A) and dp (A), satisfy the following estimates 

dn{A) < n 


dp (A) < 2 n 

We will provide upper bounds on these dimensions in terms of parameters in the model. There is a 
standard methodology to derive these estimates. We consider a volume element in the phase space, and try 
and derive conditions that will cause it to decay, as time goes forward. If A is the global attractor of the 
semigroup {S(t)} t>Q in H associated with (l)-(5) , we can define 

1 

q n (t) = sup sup - / Tr(AU(r) + F (S(r)u 0 ) o Q n (r)dr (91) 

uqEA giEH, ||< 7 i|| = l,l<i<n ^ J 0 


where 
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( 92 ) 


q n = lim sup <?„(£) 

t—yoo 

Here F is the nonlinear map in (l)-(5). also Q n is the orthogonal projection of the phase space H onto 
the subspace spanned by U\ (t), U 2 (t ), • • • , U n (t), with 


Ui(t ) = L(S(t)uo)gi,i = 1 , 2 , ..n. 


L(S(t)uo) is the Frechet derivative of the map S(t) at uo , with t fixed. Also for this model, L(S(t)uf)g = 
XJ(t) = (U, V,W, Z), where u = (u,v,w,z) is a solution to (l)-(5), <fij = (</>j...</>:•,) are an orthonormal basis 
for the subspace Q n (r)H and (C7, V, W, Z) are strong solutions to the variational equations for the reaction 
diffusion system (l)-(5) whose exact form is found in [23]. The estimate on the trace follows by standard 
means [20], and we can now state the following result as a direct application of lemma 7.3, 


Theorem 7.4 Consider the reaction diffusion equation described via, (l)-(5). The global attractor A of the 
system, for spatial dimension N < 3, is of finite dimension. Furthermore, explicit upper bounds for the 
attractors Hausdorff and fractal dimensions, are given as follows 


d H {A) < 

dp (A) < 2 


/ C(Di,a,b , c, d,a,/3) \ 2 

\ K, ) 

(1 C{Di , a, 6, c, d , a , /3) 

K~i 


|£ 1 | +1 

3 

|^| + 2 


(93) 

(94) 


7.2 Lower bound on the haursdorfF dimension of the global attractor 

We consider the reaction diffusion equation described via, (l)-(5). A spatially uniform stationary solution is 

fd 

u = w = a, v = z = —. (95) 

a 

If A is the attractor of the brusselator and L (t) is the frechet derivative of the semigroup associated to the 
brusselator, then 

dim// A > n, 

where n is the first integer satisfying 

Re\ n > 0 > Re A n _)_i, 

where {A n } represents the set of the eigenvalues of L (t). If we put 


L(t)(U 0 ,Vo,Wo,Z 0 ) = (U,V,W,Z), t > 0, 


05 I ^ 

1 

© 

> 

II 

(~(/3 + l + D 1 ) + 2uv)U +u 2 V 4- DxW, 

(96) 

dV 

m~ bAV " 

(,8 ~ 2 uv) U ~(D 2 + u 2 ) V + D 2 Z, 

(97) 

f — - 

D 3 U + (— (/3 + 1 + D 3 ) + 2 wz) W + w 2 Z, 

(98) 

f— - 

D a V +(p- 2 wz) W- (D 4 + w 2 ) Z. 

(99) 
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Replace in the linearized system, we get 




dU 

~dt 

— aAU = 

-(-P + 1 + Di)U +a 2 V + B 1 W, 



(100) 



dV 

dt 

-bAV = 

-0U - (.D 2 + a 2 ) V 

+ D 2 Z, 



(101) 



dW 

~dt 

-cAW = 

D 3 U -(-0 + 1 + D 3 )W + a 2 Z, 



(102) 



dZ 

dt 

-dAZ = 

D 4 V -0W- (£>4 ■ 

f a 2 ) Z, 



(103) 

The generator B of 

the linear semigroup L(t) 






( aA 

- (-/? + 1 +D 1 ) 

a 2 

D v 

0 

\ 


B = 



~P 

bA - (. D 2 + a 2 ) 

0 

d 2 




d 3 

0 cA 

-(-P + 1 + D 3 ) 

a 2 




V 


0 

£>4 

-0 dA — 

(Di 

+ a 2 ) ) 



Let fi\ = 0 < fi 2 < ••• < hn < ••• the eigenvalues of —A with the homogeneous Neumann boundary and 
(/?!, (f 2 ,<^nj ••• their associated eigenvectors, we try to find the eigenvectors of the operator B on the form 
$>i = (pi, qi,Ti, Si) ipi, i — 1 ,..., n, ..., with associated eigenvalues Ai, A 2 ,A n , ..., then we have 


(B - A il) = 0, 

SO 

—cipi — (—(3 -\~ 1 T T)f) — \ cf 2 D\ 0 

—(3 —bpi — {D 2 + oP) — A i 0 D 2 

D% 0 — c/i^ — (—ft + 1 + Df) — A i o? 

0 -D 4 — (3 —dfii — (D^ T of 2 ) — A i 


= 0 . 

The principal determinant of this algebraic linear system is a fourth degree polynomial of unknown A*. It 
admits a root with no positive real part if the real part of the sum of its four roots is not positive, that is 

— (a + b + c + d) Re/ii + [2 (/3 — 1 — cf 2 ) — ( D\ T D 2 + D 3 + -D 4 )] > 0 


Since 

then 


fii ~ Ki^ 


[2 (/? — 1 — a 2 ) — (Di + D 2 + Ds + Df)\ _ 2 _ 

—- J- - 7 - 77 - 1 ~ Ki* 

(cl T b T c T d) 


Thus if 


(£>i + D 2 + D 3 + D 4 ) < 2 (0 - 1 - a 2 ) 
We can state the following result 


Theorem 7.5 Consider the reaction diffusion equation described via, (l)-(5). The global attractor A of the 
system, for spatial dimension N < 3, has explicit lower bounds for its Hausdorff dimension, in particular 
there exists a universal constant K' such that 


dim// A> K' 


[ 2 ( 0 - 1 - a 2 ) - (£>1 + £> 2 + -P 3 + £>4)] 
(a T b T c T d) 


N_ 

2 
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8 Numerical Simulations 


We now carry out numerical simulations of (l)-(5). In order to explore the dynamics of the model in 2d, we 
use a finite difference method. A forward difference scheme is used for the reaction terms. For the diffusion 
terms, a standard five point explicit finite difference scheme is used. The numerical simulation is carried out 
at different time levels for two dimensional spatial model system. The system of equations is numerically 
solved over 200 x 200 mesh points , on a domain of size L x x L y , where L x = L y = 500. with spatial 
resolution Ax = Ay = 1 and time step At = ^. The initial condition used is a small perturbation about 
(2, 2.75, 2, 2.75) and the boundary conditions used are no flux Neumann conditions. Note this is fine, as the 
results hold for Neumann boundary conditions as well. 





Figure 1: Attarctor for system (l)-(5) in one dimensional domain obtained at t — 10000. 


Table 1: Parameters used in the simulations in fig 1,2,3,4 


a 

/? 

£>i 

d 2 

Ds 

£>4 

a 

b 

C 

d 

2 

5.5 

0.0126 

0.126 

0.0125 

0.125 

icr 6 

icr 6 

10“ 6 

10“ 6 


9 Attractor reconstruction 

The dynamics of the four-species Brusselator model described by (1) - (5) is visualized with the aid of 
a delay-time method. This method was proposed by Takens and Mane [19], and it ensures equivalence 
between the topological properties of the actual and reconstructed attractors. The delay method consists of 
a reconstruction performed from measurements of a single variable. Here we consider the variable u , which 
is one of the four state variables of the system. 

According to Takens and Mane [19], the dynamics of the system can be completely described by the 
time series, u(t ), which is numerically calculated at a fixed position in space. This time series can be 
embedded in an m-dimensional embedding space and reconstructed as a pseudo-trajectory through the 
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Figure 2: Spatial distribution of u,v,w,z species at t = 10000. 


following embedding vectors, which form an embedding matrix: 

Vi = (x(t 0 ), x(t 0 + t), ..., x(t 0 + (m - l)r)) T , 

U 2 = (x(to + l),x(to + l + t )),.. .,x(t 0 + l + (m- 1)t ) T , 

y s = (x(t 0 + (s - 1)1), x(t 0 + (s - 1)1 + t), ..., x(t 0 + (s - 1)/ + (to - 1)t)) t . 

(104) 

Here, r is called the “delay-time”, l is the sampling interval, and w = (m — l)r is the “window length” which 
represents the time spanned by each embedding vector. Selection of appropriate values for parameters r 
and l in the embedding procedure is particularly important for the reliability of the results. To ensure the 
equivalence between the topological properties of actual and reconstructed attractors, a formal criterion was 
proposed by Takens, namely that m > 2d + 1. This criterion relates the embedding dimension, m, to the 
attractor dimension, d. In practice, the dimension, d, of the attractor is unknown and has to be determined. 
Choosing the optimal embedding parameter, m, and the delay-time, r, is rather nontrivial. 

In this paper, we adopt the interactive technique suggested by Albano et al [1]. The authors have 
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Figure 3: The two dimensional spatial attractor of the system (l)-(5) at t = 50. 
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Figure 4: The two dimensional spatial attractor of the system (l)-(5) at t = 500. 
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Figure 5: The densities of u, v, w, z species are shown as contour plots in two dimensional plane at t = 500. 
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combined a singular-value decomposition, which leads to a set of statistically independent variables, and the 
Grassberger-Procaccia algorithm, [3], to determine the dimension of the attractors. The algorithm consists 
of four steps: 

1. Choose m and r so that the window, w, is a few times larger than the correlation time of the time 
series, u(t). A rule thumb for selecting r is to choose the time at which the autocorrelation time falls to 
about 1/e. This delay-time ensures that the embedding vectors spanning the phase space have become, 
in some sense, independent. 

2. Perform a singular-value decomposition of the embedding matrix, Y = VYU T . Singular values, cq, 
occupying the diagonal of the diagonal matrix, E, are discarded below a certain threshold, because one 
can consider them as originating from noise. In this paper, the threshold was set to be 10 -2 . 

3. Calculate the correlation integral from which the correlation dimension, d, of the attractor is then 
deduced as the limit of the slopes of the log-log plot of the correlation integrals. A test of Takens’ 
criterion is then performed; if it is not satisfied, i.e., if m < 2d+ 1, then increase m until this criterion 
is satisfied. 

4. Refine the value of the embedding dimension, m; the suitable embedding dimension is the one that 
maximizes the straight line parts in the log-log plot of the correlation integrals. 

It is well known that the addition of diffusion, to an ODE system, can lead to temporal chaos. In fact, 
the model without diffusion (the ODE case), exhibiting oscillatory dynamics for a certain parameter set, 
can be destabilized to chaos solely via adding diffusion [11]. Here, we test this observation on (1) - (5). We 
start our numerical experiments from a periodic state that exists in the absence of diffusion. The following 
parameter set yields a periodic state: 


Di = .0126, D 2 = .126, D 3 = .0125, D A = .125, a = 2, £ = 5.9. (105) 

This oscillatory dynamics is described by a limit cycle. Imposing diffusion in the four-species Brusselator 
model, that is 


a = 10“ 6 , b = 2 x 1CT 6 , c = 3 x 1CT 6 , d = 4 x 1CT 6 , (106) 

shows that the oscillatory dynamics becomes chaotic. We next reconstruct the low dimensional attractor 
from the time series of species u, which is obtained by fixing a spatial location and following the trajectory 
in time. The dimension of the reconstructed attractor, which corresponds to the slope of the correlation 
dimension of the attractor, is found to be approximately 27.54. This shows that the attractor of the full 
diffusive system is quite high. 

For the above range of parameters, we can extract a lower bound for the constant K', from Theorem 7.5. 
As 


K' 


[2(ff-l -o?)-(Di+D2 + D3+D4)] ] 2 > 2754 

(a + b + c + d) | 


one obtains after plugging in the parameters that essentially 


(107) 


K ’ loo' (108) 

Note that the fractional dimension of 27.54, indicates that the attractor is strange. To test whether or not 
the dynamics on it is chaotic, we estimated the Lyapunov exponents by the method proposed by Lai and 
Chen [9]. We find the largest exponent to be approximately equal to 2.17. This implies that the dynamics 
on the strange attractor, for this parameter set, is indeed chaotic in time. In particular, this also tells us 
that for a a parameter set where D\ ^ D 3 , D 2 ^ one may obtain a chaotic attractor, the existence of 
which becomes harder to prove, although possible, as we demonstrate. 
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